Mon. Not. R. Astron. Soc. 000, [THT2] (2008) Printed 16 December 2008 (MN I^TeX style file v2. 2) 



Properties of the ionized gas in HH 202. I: Results from 
integral field spectroscopy with PMAS* 

A. Mesa-Delgado^t, L. L6pez-Martm\ C. Esteban^, J. Garcia-Rojas^, 
and V. Luridiana^ 

• ^ Instituto de Astrofisica de Canarias, E-38200 La Laguna, Tenerife, Spain 

• ^Instituto de Astronomia, UNAM, Apdo. Postal 70-264, 04510 Mexico D.F., Mexico 

^ Instituto de Astrofisica de Andalucia (CSIC), Apdo. Correos 3004, E-18080 Granada, Spain 



Accepted X XX XX. Received XXXX XX XX; in original form XXXX XX XX 



ABSTRACT 

We present results from integral field spectroscopy with the Potsdam Multi- Aperture 
Spectrograph of the head of the Herbig-Haro object HH 202 with a spatial sampling of 
l"xl". We have obtained maps of different emission lines, physical conditions -such 
as electron temperature and density- and ionic abundances from recombination and 
collisionally excited lines. We present the first map of the Balmer temperature and of 
the temperature fluctuation parameter, t^. We have calculated the in the plane of the 
sky, which is substantially smaller than that determined along the line of sight. We have 
mapped the abundance discrepancy factor of 0^+, ADF(0^+), finding its maximum 
value at the HH 202-S position. We have explored the relations between the ADF(0^+) 
and the electron density, the Balmer and [O ill] temperatures, the ionization degree 
as well as the parameter. We do not find clear correlations between these properties 
and the results seem to support that the ADF and are independent phenomena. 
We have found a weak negative correlation between the 0^+ abundance determined 
from recombination lines and the temperature, which is the expected behaviour in an 
ionized nebula, hence it seems that there is not evidence for the presence of super-metal 
rich droplets in H ii regions. 
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1 INTRODUCTION 

The Orion Nebula is the nearest, most observed and stud- 
ied Galactic H ll region. It is an active star-formation region 
where phenomena associated with the early stages of stel- 
lar evolution such as protoplanetary disks (proplyds) and 
Herbig-Haro (HH) objects can be observed in detail. 
HH objects are bright nebulosities associated with high- 
velocity gas flows. The most prominent high-velocity fea- 
ture in the nebula is the Becklin-Neugebauer/KIeinmann- 
Low (BN/KL) complex, which contains several HH objects 
(HH 201, 205, 206, 207 and 208). But also, there are other 
important high-velocity flows that do not belong to the 
BN/KL complex, as is the case of HH 202, 203 and 204. 



* Based on observations collected at the Centro Astronomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max-Planck Institut fiir Astronomie and the Instituto de As- 
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The origin of these flows has been associated with infrared 
sources embedded within the Orion Molecular Cloud 1 South 
(see O'Del l fc Hcnncy 2008, and references therein). 
HH 202 was on e of the first HH o bjects identified in the 
Orion Nebula bv lCanto et al.l (|l980t ) in their [S ll] and [N ll] 
images as an emission line object showing two bright knots 
-the so-called HH 202-N and HH 202-S (see Figure [T])- em- 
bedded in a more extended nebulosity with a lon g concave 
form. After its discovery, the spectroscopic work bv lMeaburnI 
(I l98d ) showed that the arc-shaped nebulosity emits in [O ill] 
lines. 

Further high-spectral resolution spectroscopy in several 
emission lines ([O l], [S ll], [N ll] , [O ll], [S i l l] and Ha) 
has been performed in HH 202 bv lO'Dell etai] (|l99ll ) who 
studied HH 202-S and found two velocity components. More 
recently, we can find in the literature extensive wo rks on the 
gas kinem atics in the Orion N ebula as those by iDoi et all 
(2002 ) and lO'Dell fc Doil ||2003|). where t angential velocities 
have been measured, or jPoi et all (|2004 ). who provide with 
accurate radial velocity measurements, detecting up to three 
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Figu re 1. HST image of the central part of the Orion Nebula, which combines WFPC2 images taken in different filters llO'Dell fc Wond 
Il99d ). The white square corresponds to the field of view of PMAS IFU used, covering the head of HH 202. The separate close-up image 
on the right shows the Ha map obtained with PMAS. The original map is 16x16 pixels of l"xl" size and has been rebinncd to 160x160 
pixels. Note the remarkable similarity between the HST image and our rebinned PMAS Ha map. 



kinematic components in HH 202-S. On the other hand, the 
most detailed imaging study of the Orion Nebul a in the [S ll] , 
[O III ] and Ha lines has been carried out by lO'Dell et al.l 
l| 19971 ) with the Hubble Space Telescope (HST). These au- 
thors also detect an intense [O ill] emission in the extended 
nebulosity of HH 202 and a strong [S ll] emission at its knots. 
The extended [O lll] emission found in all the works indicates 
that the main excitation mechanism of HH 202 is photoion- 
ization, most probably from the brightest and hottest star 
of the Trapezium cluster, 0^ Ori C. This particular property 
-not very common in other HH objects, which are usually 
ionized by shocks- permits to derive the physical conditions 
and chemical abundances of the high velocity gas associated 
with HH 202 making use of the standard techniques for the 
study of ionized nebulae. 

Our group has studied in detail the chemical composition of 
the Orion Nebula in some previous pa pers l|Esteban et al.l 
1 19981 . 12004 iMesa-Delgado et al.1 |2008| ) and is especially 
interested in studying the so-called abundance discrep- 
ancy (AD) problem in H ll regions, which is the dis- 
agreement between the abundances of the same ion de- 
rived from colhsionally excited lines (CELs) and recom- 
bination lines (RLs). From intermediate and high resolu- 
tion spectroscopy of Galactic and extragalactic H ll re- 
gions, our group has found that the 0'^'''/H^ ratio calcu- 
lated from RLs is always between 0.1 and 0.3 dex higher 
than the value obtained from CELs for the same region (see 
iGarcfa-Roias fc Estebanll2007l ). The results obtained for H ll 
regions are quite different to those obtained for planetary 
nebulae (PNe), where the AD shows a much wider range of 
values and could be substantially larger in some objects {e.g. 
iLiu et aLllioOol. l2006l: pTsainis et all 120031 . [2004). As argued 
in lGarci'a-Roias fc EstebanI (|2007l '). the AD problem in H ll 
regions seems to be consistent with the predicti ons of the 
tempe rature fluctuaction paradigm proposed by IPeimber^ 
l|l967t) and parametrized by the mean square of the spa- 
tial distributions of temperature, the so-called param- 
eter. In this scenario, the AD would be produced by the 
very different temperature dependence of the emissivities 



of both kinds of lines. In any case, the existence and ori- 
gin of such temperature flu ctua tions are still controv ersial. 
iTsamis fc PequignotI (120051 ) and lStasihska et all (|2007l ) have 
proposed a different hypothesis in order to explain the AD 
in H II regions. Base d on t he model for heavy-element mix- 
ing of Tcnorio-Tagy (|l996l ). the presence of cold metal-rich 
droplets of supernova ejecta still not mixed with the ambi- 
ent gas of the H ll regions would produce most of the RL 
emission while GEL emission would be produced by the am- 
bient gas, wich would have the typical electron temperature 
of an H II region and the expected composition of the lo- 
cal interstellar medium. In such case, the abundances from 
GELs would be the truly representative ones of the nebula. 
Several authors have studied the spatia l distribution of the 
physical conditions in the Orion Nebula. [OT>ell et al.1 (|2003| ) 
have obtained a high-spatial resolution map of the elec- 
tron temperature in a field centered at the southwest of 
the Trapezium cluster from narrowband images taken with 
the HST. They found small spatial scale temperature vari- 
ations, which seem to be in quantitative agreement with 
th e parameter determined using the AD factor found 
bv lEsteban et al.1 (|2004D . On the other hand, iRubin et al.l 
|2003l ) have taken long-slit spectroscopy making use of ST IS 
at the HST in several positions and do not find substan- 
tial spatial variations or gradients of the physical condi- 
tions a long their slits positions. More recentlv. lSanchez et al.l 
l|2007l ) have obtained a mosaic of the Orion Nebula from in- 
tegral field spectroscopy with a spatial resolution of 2'.'7, 
but lacking a confident absolute flux calibration due to bad 
transparency during the observations. Their electron tem- 
perature map shows clear spatial variations and the electron 
density map is very rich in substructures, some of which 
could be related to HH objects. A more detailed view - 
but based on long-s li t spec troscopy- has been obtained by 
|m esa-Delgado et al.l l|2008h . who obtained spatial distribu- 
tions -at spatial scales of 1'.'2- of electron temperature, den- 
sity, line intensities, ionic abundances and the AD factor 
(ADF, defined as the difference between the abundances de- 
rived from RLs and GELs) for five slit positions crossing 
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different morphological zones of the Orion Nebula, finding 
spikes in the distribution of the electron temperature and 
density which are related to the position of proplyds and 
HH objects. In particular, these authors found that the ADF 
shows larger values at the HH objects. 

The use of the integral field spectroscopy is still in its in- 
fancy in the study of ionized nebulae and the number of 
works available is st ill rather small. In the case of PNe, 
iTsamis et al.l (|2008bl) have carried out the first deep study 
using this technique with FLAMES at the VLT taking spec- 
tra of three PNe of the Galactic disc and covering the spec- 
tral range fr o m 396 4 to 5078 A. For the Orion Nebula, 
iTsamis et al.l l|2008ah have presented preliminary results of 
integral field spectrocopy with FLAMES for several proplyds 
with the aim of studying the behaviour of the physical con- 
ditions, chemical abundances and the ADF in and outside 
these kinds of objects. The first results show that the tem- 
perature measured in the proplyds can be affected by col- 
lisional deexcitation as has been also suggested in p revious 
works (|Rubin et all |2003| : iMesa-Delgado et al.l boosi ) . 
The main goal of the present paper is to use integral field 
spectroscopy at spatial scales of 1" in order to explore the 
effect of HH objects on the derived ADF considering the re- 
suhs of iMesa-Dekado et all (j2008l 'l. 

In ^ we describe the observations obtained with PMAS 
and the reduction procedure. In ij3] we describe the emis- 
sion line measurements and the reddening correction as well 
as some representative bidimensional maps of those quan- 
tities. In fj4] we describe the determination of the physical 
conditions, chemical abundances -from both kinds of lines, 
CELs and RLs- and the ADF of O^^ as well as their cor- 
responding maps. In ^we calculate the parameter from 
different methods and present the first bidimensional map of 
this quantity. Also, we show the possible correlations among 
the ADF and the rest of properties determined. Finally, in 
3n]we summarize our main conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

HH 202 was observed on 2007 October 14 at Calar Aho 
Observatory (Almeria, Spain), using the 3.5m Telescope 
with the Potsda m Multi- Aperture Spectrometer (PMAS, 
iRoth et al.ll2005l ) . The standard lens array integral field unit 
(IFU) of 16"xl6" field of view (FOV) was used with a sam- 
phng of 1". Most of the optical range was covered with the 
V600 grating using two grating rotator angles: —72, cover- 
ing from 3500 to 5100 A; and —68, covering from 5700 to 
7200 A. The effective spectral resolution was 3.6 A. The po- 
sition of the IFU covering HH 202 is shown in Figure[T]where 
we can see the two bright knots HH 202-N and HH 202-S 
together with an Ha image obtained with PMAS rebinned 
to 160x160 pixels using the rebin function of idl which 
performs an expansion of the original map through a linear 
interpolation. The blue and red spectra have a total integra- 
tion time of 1800 and 1000 seconds, respectively. Additional 
short exposures of 10 seconds were taken in order to avoid 
saturation of the brightest emission lines. Calibration images 
were obtained during the night: arc lamps for the wavelength 
calibration and a continuum lamp needed to extract the 256 
individual spectra on the CCD. Observations of the spec- 
trophotometric standard stars BD +28°4211, Feige 110 and 
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Figure 2. Section of a PMAS spectra around O II lines of multi- 
plot 1 corresponding to the spaxel position (-5.5, -6.5). 

G 191-B2B (|Okelll990l ') were used for flux cahbration. The 
error of this calibration is of the order of 5%. The night was 
photometric and the typical seeing during the observations 
was 1". 

The data were reduced using the iraf reduction package 
SPECRED. After bias subtraction, spectra were traced on the 
continuum lamp exposure obtained before each science ex- 
posure, and wavelength calibrated using a HgNe arc lamp. 
The continuum lamp and sky flats were used to determine 
the response of the instrument for each flber and wavelength. 
Finally, for the standard stars we have co-added the spectra 
of the central flbers and compared them with the tabulated 
one-dimensional spectra. 

We have noticed the effect of the differential atmospheric 
refraction (DAR) in the monochromatic images of HH 202 
obtained for Balmer lines at different wavelengths reaching 
the value of ~l'.'3 between Ha and Hll. We have measured 
offsets between all Balmer line images, and shifted with re- 
spect to Ha. The maximum coincident FOV resulting in the 
whole wavelength range is 15"xl5". All the maps involving 
emission line ratios analysed in this paper have been cor- 
rected for DAR. 

Due to the relatively low spectral resolution of the obser- 
vations, Hg 4358 A telluric emission was somewhat blended 
with [O III] 4363 A and a proper sky subtraction was neces- 
sary. We have taken a median value of the Hg 4358 A line 
flux (9.22 xlO"^'' erg cm~^ s~^), which is in some regions 
approximately two orders of magnitude lower than the au- 
roral line flux. 

In previous works the background nebular emission in 
Orion proplyds has been de termined, either using echelle 
data (|Hennev &: 0'DellllT999l ) or integral fleld spectroscopy 
(IVasconcelos et al. ||2005'). In our case, the low spectral res- 
olution of the data does not allow to separate the HH 202 
emission from the background. Furthermore, as pointed out 
in both cited papers, this task is very challenging due to 
small-scale inhomogeneities in the nebular emission. As a 
result, only rough estimates can be given. Using the method 
proposed bv iVasconcelos et al . (200^) , we have assumed that 
the background emission is represented by the emission in 
the northeast corner of the FOV. The resulting background 
contribution is thus estimated to lie between 40% and 50% 
in knot HH 202-S, and between 55% and 65% in the rest 
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of HH 202. If such a correction were applied to our data, 
we would obtain much larger densities, lower temperatures 
from the [N ll] line ratio and higher temperatures from the 
[O III] line ratio with respect to the results presented in next 
sections. These results do not change the final conclusions 
of the paper because we are interested in how the presence 
of an HH object affects to the ADF when low spectral reso- 
lution observations are made and the contributions from the 
gas flows and the nebular background are mixed. In Paper II, 
we perform a complete analysis of the knot HH 202-S based 
on echelle spectrophotometry. The high resolution of these 
data has allowed us to resolve the kinematic component as- 
sociated with the gas flow and the one associated with the 
background emission (in a small area included in the PMAS 
field). This analysis indicates that the background emission 
ranges from 30% in the case of [N ll] 5755 A to 60% for [O ill] 
5007 A. In the case of H I lines, the background emission is 
about 50%. 



3 LINE MEASUREMENTS AND REDDENING 
CORRECTION 

The emission lines considered in our analysis are the fol- 
lowing: a) Balmer lines, from Ha to Hll, which are used 



to compute the reddening correction and correct the DAR 
shift (H8 and HIO were not used because they suffer from 
line-blending); b) CELs of various species, which are used to 
compute physical conditions and ionic abundances; c) Faint 
RLs of C II and O ll which are used to derive the ionic abun- 
dances and to compute the ADF(0^^). In Figure [51 we show 
an example of the spectra of a given spaxel where we can 
see the blend of O ll lines at 4649 and 4651 A. 
Line fluxes were measured applying a single or a multi- 
ple Gaussian profile fit procedure between two given lim- 
its and over the local continuum. All these measurements 
were made with the SPLOT routine of iraf and using our 
own scripts to automatize the process. The errors associ- 
ated w ith the line fiux measurem ents were determined fol- 
lowing [Mes^DeigadoelaD l|2008l '). The final error of a line 
was computed as the quadratic sum of the error in its flux 
measurement and the error in flux calibration. In order to 
avoid spurious weak line measurements, we imposed three 
criteria to discriminate between real features and noise: 1) 
Line intensity peak over 2.5 times the sigma of the contin- 
uum; 2) FWHM(H i)/1.5 < FWHM(A) < 1.5xFWHM(H i); 
and 3) F(A) > 0.0001 x F(H/3). 

All line fluxes of a given spectrum have been normalized to 
H/3 and Hq for the blue and red range, respectively. To pro- 
duce a final homogeneous set of line ratios, all of them were 
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Figure 4. HH 202 PMAS maps with Ha contours overplotted. a) c(H/3), b) Ha heliocentric velocity. 
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Figure 5. Physical conditions maps with Ha contours overplotted: a) ne{[S ll]), b) Te{[0 III]), c) Te([N ll])/Te{[0 III]) ratio and d) 
T41N II]). 
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re-scaled to H/3. The re-scaling factor used in the red spectra 
was the theoretical Ha/H/3 ratio for the physical conditions 
of Te= 10000 K and ne= 1000 cm"^ 

In Figure H we show [O ill] 5007 A, O ii 4650 A, [N ii] 
6584 A and [Fe ill] 4881 A emission line maps with Ha con- 
tours overplotted. In Figures |3K andlSh, it is possible to see 
that the [O lll] 5007 and O ll 4650 A spatial distributions 
are quite similar. In particular, the O ll flux map seems to 
peak at the position of the apex of HH 202, the zone also 
known as HH 202-S -spaxels (-5.5,-6.5) and (-6.5,-4.5)-. A 
similar spatial distribution of both kinds of lines was al- 
re ady observed using lo ng-slit spectra in the Orion Nebula 
bv lMesa-Delgado et aP (2008), as well as in the three PNe 
observed bv lTsamis et al.. (.2008b ). Figure |3] also shows that 
the spatial distributions of [N ll] and [O ill] are rather differ- 
ent. The [N ii] map shows an enhancement in the bow shock 
of HH 202. HH objects are characterized by their strong 
emission in [Fe ill] lines. In Figure [3ji we can see the spatial 
distribution of [Fe ill] 4881 A where the maximum emission 
is strongly concentrated at HH 202-S, which coincides with 
the maxima in radial velocity and electron density as we 
can see in Figures |4]3 and[5K. This high [Fe ill] emission at 
HH 202-S may be related to shock destruction of dust grains 
or to the particular excitation conditions of the gas in the 
knot. 

The reddening coefficient, c(H/3), has been obtained by fit- 
ting the observed H7/H/3, H(5/H/3, H9/ H/3 and HI 1 /H/3 ra- 
tios to the theoretical ones predicted bv lStorev fc Hummer] 
l|l995l) for ne = 1000 cm"^ and Te = 10000 K. We have 
used the red dening function , ,f (A) , normalized to H/9 de- 
termined by 'Blagravc et al.l l|2007l ) for the Orion Nebula. 
T he use of this extinction la w instead of the classical one 
of lCostero fc PeimbertI l|l970l ) produces c(H/3) values about 
0.1 dex higher and dereddened fluxes with respect to H/3 
about 3% lower for lines in the range 5000 to 7500 A, 4% 
higher for wavelengths below than 5000 A and 2% for wave- 
lengths above 7500 A. The final adopted c(H/3) value for 
each spaxel is an average of the individual values derived 
from each Balmer line ratio weighted by their correspond- 
ing uncertainties. The typical error of c(H/3) is about 0.14 
dex for each spaxel. The resulting extinction map is shown 
in Figure 2^. The extinction coefficient varies approximately 
from 0.4 to 0.7 dex and reaches the highest values -between 
0.8 and 0.9 dex- in a low surface brightness area at the north 
of HH 202. Excluding this zone -which shows the largest line 
intensity uncertainties- we obtain a mean value 0.56±0.14 
dex. 

We have compared our c ( H/3) v alues with those obtained 
bv lO'DeU fc Yusef-Zadehl l|2000l ) and iMesa-Delgado et all 
(|2008t ) for the area of HH 202. The first authors have ob- 
tained consistent c(H/3) values from the Hq/H/? line ratio 
using calibrated HST images and from radio to optical sur- 
face brightness ratio. HH 202 is locat ed between the contour 
lines of 0.2 and 0.4 dex on t he maps oflO'Dell fc Yuscf-Zadeh' 
(2000). On the other hand. lMesa-Delgado et" (^008), from 
long-slit spectroscopy passing precisely through HH 202-S, 
derive a c(H/3) value of about 0.4±0.1 dex. Therefore, we 
find that our mean value of the reddening coefficient is con- 
sistent with the previous determinations taking into account 
the errors and the different extinction laws used. 
The heliocentric velocity map obtained from the centroid of 
the Ha line profile is shown in Figure UJs. We observe that 
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Figure 6. Section of a PMAS spectrum of the brightest spaxel 
at the apex (-5.5,-6.5). The fitted continuum contribution to the 
right side of the BJ was substracted to the whole range showed. 



the gas around the head of HH 202 shows a velocity of about 
—5 km s~^ and that the most negatives values are reached 
at HH 202-S position, where velocities of —35 km s"'^ are 
measured. These values of the velocity of HH 202-S agree 
within the errors with that determined bv lDoi et al.l (|2004l ) 
of — 39±2 km s~^ for this position. However, due to the lower 
spectral resolution of our data, the velocities we measure 
correspond to the blend of the emission of the background 
nebular gas and the gas flow of the HH object. 



4 RESULTS 

4.1 Physical conditions 

Nebular electron densities, n^, and temperatures, Te, 
have been derived from the usual GEL ratios -[S 11] 
6717/6731 for ne, and [O iii] (4959-^5007) /4363 and [N 11] 
(6548-|-6584)/5755 for Te, a nd using the temden task of the 
NEBULAR package of IRAF (Shaw & Dufour 1995) with up- 
dated atomic data (see iGarcia-Rojas et al.,,,20p5,). F ollowing 



the same procedure as iMesa-Delgado et al.l (120081) for the 



determination of the physical conditions, we have assumed 
an initial Te= 10000 K to derive a first approximation of 
rie, then we calculate Te{[0 ill]) and re([N 11]), and iterate 
until convergence. We have not corrected the observed in- 
tensity of [N 11] 5755 A for contribution by recombination 
when determining Te([N 11]) because this contribution is ex- 
pected to be rather small in H 11 region s with high ionizatio n 
degree such as the Orion Nebula (e.g. lEsteban et al]l2004l '). 
Bidimensional maps of We, Te{[0 ill]), Te([N ll])/T'e([0 ill]) 
ratio and re([N 11]) are shown in Figures [5^, [S]d,[SJ;, and[S}l, 
respectively. 

In Figure [5^ we can see that the highest densities are just 
reached at the positions of HH 202-N and HH 202-S with 
val ues of about 10000 cm""^. Similar values were obtained 
by iMesa-Delgado et all (120081 ) for HH 202-S. For the gas 
outside the HH object, we find a We of about 4000 cm"'' 
and the lowest densities (~2500 cm~^) at the low surface 
brightness region at the northwest corner of the FOV. The 
typical errors in our density determination are between 200 
and 250 cm"^. 

Figure [5]d shows that Te([0 ill]) is rather constant in the 
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Figure 7. Balmer temperature map with Ha contours overplot- 
ted. The black rectangle on the northwest corner corresponds to 
an area masked due to the bad determination of Te(Bac). 



whole FOV. However, the highest values of re([0 ill]) trace 
the form of HH 202 -with an almost constant temperature 
of 8500 K- reaching the peak value at the HH 202-S posi- 
tion where shock heating should be maximum. On the other 
hand, re([N ll]) shows a much wider range of values across 
the FOV, but its spatial distribution is completely differ- 
ent to that of Te{[0 III]). The maximum values of re([N ll]) 
are found just behind HH 202-S or ahead of the northwest 
border of the bow shock of HH 202-N (see Figure [5jl). The 
lowest values are found in the inner part of HH 202 (~9500 
K). In Figure [St, the re([N l\])/Te([0 ill]) ratio illustrates 
that Te([N ll]) is always higher than Te([0 ill]) but the ratio 
tends to be closer to one in the HH object and especially at 
HH 202-N where the density is higher. The uncertainties of 
the temperatures are about 600 K for Te([N ll]) and 250 K 
for Te{[0 III]). 

The Balmer continuum temperature, Te(Bac), of the ion- 
ized gas has also been determined for most of the individual 
spaxels of the FOV; this is the first bidimensional Te(Bac) 
map ever obtained for an ionized nebula. This temperature 
depends on the ratio between the Balmer jump (BJ) and the 
intensity of a given Balmer line. We have measured the BJ 
subtracting the continuum fitted at both sides of the Balmer 
discontinuity at 3646 A (see Figure |6]) . The continuum fit- 
ted at both sides of the BJ is defined over spectral ranges 
free of emission lines covering the possible maximum base- 
line. The errors associated with the continuum subtraction 
are included in the discontinuity measurements. Finally, we 
have computed Te(Bac) in K from the r atio of the BJ to the 
Hll flux using the relation proposed bv lLiu et all l|200lh : 

T,(Bac) = 368 x (l-f 0.259y+ 3.409y^+) ~ ' (1) 

where and 2/"*"+ correspond to the He"'"/H"'" and He^"'"/H"'" 
ratios, r espectively. We have assumed for the value de- 
rived by lEsteban et all l|2004l ) for their slit position at the 
center of the Orion Nebula and y"*"^ = due to the lack 
of He II lines in our spectra. The effect of the interstellar 
extinction was taken into account in the calculations. 



The re(Bac) map is shown in Figure [7] The Te(Bac) error 
was calculated by error propagation in the equation (1), and 
it amounts to 22% on average (i.e., 1100 to 1900 K). We can 
see that the lowest temperatures (~5500 K) are reached at 
the east half of the FOV, while the west half shows values of 
about 7000 K, encompassing HH 202-S and HH 202-N. From 
the comparison of Figure [7] and the rest of the maps, we can 
see that the spatial distribution of Te(Bac) is something sim- 
ilar that of ^^([O III]): both indicators show the maximum 
values at HH 202-S, HH 202-N and the arc connecting both 
features. 

4.2 Chemical Abundances 

The IRAF package nebular has been used to derive ionic 
abundances of N"*", O"^ and O^"^ from the intensity of CELs. 
We have assumed no temperature fluctuations in the ionized 
gas (t^= 0) and a two-zone scheme: re([N ll]) for N^ and O"*" 
-the low ionization potential ions- and Te([0 ill]) for O^"*". 
The errors in the ionic abundances have been calculated as a 
quadratic sum of the independent contributions of errors in 
flux. Tie and Te. The spatial distribution of the O^""" and O""" 
abundances are presented in Figures[SjD and[S]i, respectively. 
In these flgures, we can see that the spatial distributions of 
both abundances are completely different. The O^/H"*" ratio 
reaches the highest values just on the bow shock, whereas 
0^+/H+ shows an inverse behaviour, reflecting the different 
ionization structure of the bow shock and the bulk of the 
background nebular gas. The higher 0'''/0'^^ ratio at the 
bow shock is probably produced by the higher densities in 
the shock gas, which increase the recombination rate of O'^"'" . 
The abundace map of N^ -not presented here- shows a very 
similar distribution to the O^ map. 

The high signal-to-noise ratio of our spectra has permitted 
to detect and measure the weak RLs of O ll and C ll in all 
the spaxels of the PMAS FOV. These RLs have the advan- 
tage that their relative intensity with respect to a Balmer 
line is almost independent on Te and n^, largely unaffected 
by the possible presence of temperature fluctuations. In Fig- 
ure [2] we show a section of the spectra at spaxel (-5.5, -6.5) 
around the O 1 1 lines. 

We have measured the blend of the O ll 4649 and 4651 A 
lines -the brightest individual lines of the O ll multiplet 1- in 
most of the spaxels of the FOV. In order to determine the 
O^""" abundance fr om this blend, we hav e com bined equa- 
tions (1) and (4) of lPeimbert fc Peimberl l|2005h estimating 
the expected flux of all the lines of O ll multiplet 1 flux with 
respect to H/3, /(Ml O ll)//(H/3). These relations were used 
in order to correct for the NLTE effects in the relative inten- 
sity of the individual lines of the multiplet, though these are 
rather small in the Orion Nebula due to its relatively large 
Ue. Then, the O^^/H""" ratio is calculated by: 

O^ ^ Am^ ae//(J//?) J(M1 O II) 

H+ 4861 Qe//(M1) I{HI3) ' ^ ' 

where a^ff{HI3) and ae//(Ml) are the effective recombi- 
nation coefficients for H/3 and for the O ll multiplet 1, re- 
spectively, and Ami = 4651.5 A the representative mean 
wavelength of the whole multiplet. We can also obtain the 
C^^ abundance using an equation analogous to (2) but us- 
ing the specific quantities for this ion. Both abundances have 
been calculated assuming T'e([0 ill]) and with the effective 
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Figure 8. Ionic abundance maps with Ho contours overplotted: a) 12 + log(0"''^/H+) from RLs —the black rectangle on the northwest 
corner corresponds to an area masked due to the bad determination of the abundance, b) 12 + log(0"'"+/H+) from CELs, c) 12 + 
Iog(C++/H+) from RLs and d) 12 + log(0+/H+) from CELs. 



recom bination coefficients available in the literature ([Storevl 
1 19941 for assuming LS coupling and iDavev et al.ll2000l 
for C^+). Typical ionic abundance errors are 0.10 dex for 
and 0.08 dex for C^+. The ionic abundances obtained 
from RLs are shown in Figures [8^ and|8j;. 
The spatial distributions of the O^"'' abundance obtained 
from CELs and RLs show that the lower values are found 
at the positions of HH 202-S and HH 202-N and the arc 
connecting both features. In contrast, the O^^ abundance 
map obtained from CELs shows especially large values in 
the zone inside the HH object at the bottom of the FOV. 
For a given spaxel, the O^^/H"*" ratios obtained from RLs are 
always higher than those obtained from CELs, with a typi- 
cal difference of about 0.10 dex. Finally, the C^^ abundance 
map shows rather constant values of about 8.35 dex, with 
its minimum value also found at HH 202-S. The qualitative 
behaviour of this map is rather similar to the O^^ abun- 
dance spatial distribution. 

For each spaxel, the ADF(0^^) has been calculated from 
the difference between the O^^ abundances obtained from 
RLs and CELs. In Figure [9] we show the spatial distribution 
of the ADF(0^'*'), which shows the highest values at and 
around HH 202-S (up to 0.23 dex), the rest of the spaxels 
show values between 0.05 and 0.20 dex with a typical error 
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Figure 9. ADF(02+) map with Ha contours overplotted. The 
black rectangle on the northwest corner corresponds to an area 
masked due to the bad determination of the O^"*" abundance from 
RLs. 
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value of 0.12 dex. Excluding the masked area showed in Fig- 
ure |9l we have calculated an ave rage value of 0.13 dex with a 
standard deviation of 0.05 dex. iMesa-Delgado et al.l (|2008l ) 
also found a peak of the ADF(0^^) of about 0.30 dex at the 
same zone of HH 202. Both values of the discrepancy are 
consistent within the uncertainties, but the value obtained 
by those authors is, in principle, more reliable because of 
the better signal-to-noise ratio of their spectra. 

5 DISCUSSION 

5.1 Temperature fluctuations 

One of the possible explanations for the AD problem in 
ionized nebulae i s provided by the presence of t emperature 
fluctuations (e.g . iGarci'a-Roias fc Estebanll2007l ). Following 
|Peimber3 (|l967l ). the temperature fluctuation over the ob- 
served volume of a nebula can be parameterized by 

To = ^4 ^ ' 3 

where To is the average temperature of the nebula and rii is 
the ion density, and 

TiJn.mdV ' ^' 

where is the mean-square electron temperature fluctua- 
tion. 

There are sever al methods to obtain and To along the 
line of sight (see lPeimbert et"aL I bOO^ ). One possibility is to 
compa re the Te obtain ed from independent methods. Fol- 
lowing iPeimberO l|l967f ) , the temperature derived from the 
BJ depends on and To as: 

TeiBac) = ro(l - l.70t^) , (5) 

and the te mperatures derived from th e ratio of CELs can be 
written as (jPeimbert fc Costerd[l969h 

T,(h) = ro(h)[n-i(|^-3)t^(h)] , (6) 

and 

Te(l)=To(l)[l + i(^-3K(l)], (7) 

where re(l) and Te(h) are the electron temperatures for the 
low -re([N ll])- and high -re([0 ni])- ionization zone, re- 
spectively. On one hand, we can obtain the t^{0'^^) and 
To(0'^^) combining equations (5) and (6). On the other 
hand, i2(N+) andTo(N+) can be obtained combining expres- 
sions (5) and (7). To obtain the mean square temperature 
fluctuation for the entire volume, , we have weigthed the 
relative importance of f^(N^) a nd ^^(O^^) in the ob served 
volume by using equation 16 of IPeimbert et al] (|2002i ). as- 
suming that in this region and O^ coexist in the same 
volume. In Figures [Tok. \TOk and [lOt;, the t^(0^+), t^(N+) 
and < > maps are presented. Their associated error maps 
are also shown in Figures [TUb . llUt l and II Of . 
The inspection of Figures [TOb . llOfc and llOb indicates that the 
spatial distribution of values are rather similar, showing 
somewhat lower values at the positions of HH 202-S and 
HH 202-N. We have calculated the mean value of < > in 
the area where this quantity has been determined weighting 



the individual values of the fluctuation for each spaxel by 
its associated error, obtaining < t'^ > ^ 0.061±0.022. 
Another possibility to obtain the parameter is assuming 
that the ADF(0^^) is produced by the presence of temper- 
ature fluctuations. Following this as sumption, we ha v e used 
equations (9), (10) and (11) from iPeimbert et all (|2004 ) 
to derive the associated 4^(0^"'") for each spaxel that -as 
expected- shows a spatial distribution analogous to the ADF 
map shown in Figure In this case, the mean value of is 
only slightly lower than that obtained from the data of Fig- 
ure llOb . A remarkable result is that the maps obtained 
from the two different methods -comparison of tempera- 
tures and from the ADF(0^^)- are qualitatively different 
suggesting, in this case, that there is not a genetic relation- 
ship between and ADF(0^+). 

Finally, we have determined the in the plane of the sky, 
t^, us ing equations (11) and (12) of Mcsa-D elgado et al.l 
l|2008l ). averaging the bidimensional point-to-point varia- 
tions of re(0^+) and Te(N+). The values obtained are 
tA(0^+) - 0.0004 and f4(N+) 0.0023, in agree ment with 
the long-slit results bv lMesa-D elgado e t al.l (|2008t ). and sub- 
stantially lower than the values obtained with the meth- 
ods described in the previous paragraphs, which correspond 
to the ffuctuations along the line of sight. 

5.2 Correlations among the ADF and other 
nebular properties 

In this section, we investigate some relationships among 
physical properties derived in the HH 202 PMAS fleld. In 
Figure [TT] the possible correlations among the ADF(0^^), 
c(H^), and ionic abundances determined from 
CELs and RLs are presented. 

The dependence of the ADF(0^"'') with respect to c(H/3), 
Tie, Te([0 III]), Te(BJ), t^{0^+) aud 0^+/0+ ratio is shown 
from Figures [TTb toOTt. Figure [TTb -ADF(0^+) vs. c(H/3)- 
does not show a clear trend, indicating that the ADF does 
not depend on the amount of dust prese nt in the line of sight 
of each spaxel of the PMAS fleld. As 'Mesa-Delgad o et al.l 
(J008) also obtained from their long-slit data, there seems 
to be no correlation between the ADF and Ue (see Figure 
lllb ) as well as between the ADF and the ionization degree 
of the gas (see Figure . 

In Figures lllb and IllH . we can see that positive correla- 
tions are obtained between ADF(0^"'') vs. Te([0 ill]) and 
ADF(0^+) vs. Te(BJ), respectively. This is expected be- 
cause these three quantities show larger values at HH 202-S. 
As we can see in Figures [TTb and lllt l. these correlations do 
not seem very reliable considering the very narrow baseline 
of Te{[0 III]) and Te(BJ) covered by the data and the large 
error bars with respect to the baseline. 

A rather unclear negative correlation between the 
ADF(0^^) and the fluctuation parameter t^{0^'^) along the 
line of sight can be guessed in Figure [TTfe. However, the large 
error bars in both quantities do not permit to establish any 
correlation. This result suggests the ADF and may be in- 
dependent. 

In Figures lllk . Illh and IIH . we represent the O^^ abun- 
dances from RLs, O'^"'' abundances from CELs and C^"^ 
abundances from RLs vs. Te{[0 m]), respectively. In the case 
of the abundances from RLs, we can observe negative cor- 
relations with rather low Spearman correlation coefficients 
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Figure 10. parameters and their respective error maps in percentage witli Ha contours overplotted: a) and b) t^(0^^) and its error, 
c) and d) t-^(N+) and its error, e) and f) and its error. The white areas are where the error is higher than the value. The black 
rectangle on the northwest corner corresponds to an area masked due to the bad determination of Te(Bac). 
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(p = —0.42 and p = —0.28 for Figures [TTk and lilt , respec- 
tively). This apparent trend is mainly produced by the be- 
haviour observed at HH 202- S and HH 202-N, which are 
precisely the zones where the O^^ and C^^ abundances ob- 
tained from RLs are lowe r and Te([0 i l l]) is h igher. From 
integral field spectroscopy, iTsamis et all (j2008bl ) have found 
rather clear positive correlations between ionic abundances 



obtained from RLs and Te in three PNe, exactly the op- 
posite trend we observe in this FOV of the Orion Nebula. 
These authors argue that their results provide evidence for 
the existence of two distinct components of highly ionized 
gas at very different temperatures. However, in this case, we 
observe the expected natural behaviour of an ionized nebula 
-the O'^"'' abundance from RLs increases when the temper- 
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Figure 11. Correlations among different nebular properties obtained from PMAS maps. Prom a) to f) the ADF{0^+) is represented 
with respect to c(H/3), Ue, Te([0 m]), Te(BJ), (^[O ill] and log(02+/0+). The mean error bar in the ADF(02+) is only shown in a). 
Prom g) to i) the logarithmic abundances of from RLs, from CELs and from RLs are represented vs. Te([0 III]). 



ature decreases- due to the fact that O^^ is a dominant 
coolant in the gas phase and the presence of spatial changes 
of the O^'^ /O^ ratio across the FOV. Therefore, the pres- 
ence of metal-rich droplets does not seem to be supported 
by our results. On the other hand, the ratio from 

CELs shows a much tighter correlation with Te([0 ill]) due 
to the lower uncertainties of these abundances. We have fit- 
ted the following linear relation: 

02+ 

12 + log(-fj-) = (11.4 ± 0.1) + (-3.6 ± 0.1)T4 (8) 

where T4 is Te([0 ill]) in units of 10* K and we have obtained 
a linear correlation coefficient r = —0.90. 

6 CONCLUSIONS 

In this paper, we present results from integral field spec- 
troscopy of an area of 15" x 15" covering the head of HH 202 
in the Orion Nebula. The FOV comprises the bright regions 
known as HH 202-S and HH 202-N. We have obtained maps 
of relevant emission line ratios, physical conditions and ionic 
abundances, including O^"*" and C^"*" abundances determined 



from recombination lines (RLs). Additionally, we have ob- 
tained maps of other interesting nebular properties, as the 
temperature fluctuation parameter, t^, and the abundance 
discrepancy factor of O^"*", ADF(0^^), which is defined as 
the difference between the 0^^/H+ ratios determined from 
RLs and coUisionally excited lines (CELs). 
We find that the flux distributions of the [O in] and O 11 
lines are rather similar and that the HH object is compara- 
tively much brighter in the lines of low-ionization potential 
ionic species and especially in [Fe ill]. The Ue is about 4000 
cm~^ in most of the FOV and higher -about 10000 cm""^- 
at HH 202-S and HH 202-N. The Te([0 ill]) map shows 
a narrow range of variation, but the values are higher at 
HH 202-S. On the other hand, re([N 11]) shows larger varia- 
tions and a very different spatial distribution, being higher 
at the northern and eastern edges of HH 202-N and HH 202- 
S, respectively, a likely consequence of the ionization strati- 
fication or the presence of some shock excitation in the knots 
of the HH object. We have obtained -for the first time in an 
ionized nebula- the Te(Bac) map, which follows closely that 

of Te([0 III]). 

The 0"'"/H"'" ratio map reaches the highest values just on the 
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arc that delineates the north of HH 202 from the east and 
west edges of the FOV, whereas the map of O^"'' abundance 
obtained from CELs shows an inverse behaviour, probably 
produced by the higher densities in the shock gas that in- 
crease the recombination rate of O'^"''. The spatial distribu- 
tions of the O^"'' abundance obtained from CELs and RLs 
agree in showing the lower values at the positions of HH 202- 
S and HH 202-N and the arc connecting both features. How- 
ever, for a given spaxel, the 0^"*'/H^ ratios obtained from 
RLs are always about 0.10 dex higher than those obtained 
from CELs. The map of the ADF(0^+) shows the highest 
values at and around HH 202-S. 

We have determined -for the first time in an ionized nebula- 
the map of the FOV from the comparison of re(Bac) and 
the other Te values determined from coUisionally excited line 
ratios, finding that it does not match the ADF(0'^''') map. 
This result supports that the AD and temperature fiuctua- 
tions are independent phenomena. 

Finally, we have found weak correlations between the 
ADF(0^^) and Te values of the different spaxels of the 
FOV and between the 0'^'''/H+ ratios obtained from RLs 
and CELs and T'e([0 ni]). The negative correlation between 
the O'^"'' abundances obtained from RLs and T'e([0 ill]) does 
not suppo rt the predictions of t he che mica l inhomogeneity 
models of iTsamis fc PeauignotI l|2005h and IStasihska et al.l 
(|2007l ). 
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